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Abstract — We propose a selection region based multi-hop 
routing protocol for random mobile ad hoc networks, where 
the selection region is defined by two parameters: a reference 
distance and a selection angle. At each hop, a relay is chosen 
as the nearest node to the transmitter that is located within 
the selection region. By assuming that the relay nodes are 
randomly placed, we derive an upper bound for the optimum 
reference distance to maximize the expected density of progress 
and investigate the relationship between the optimum selection 
angle and the optimum reference distance. We also note that 
the optimized expected density of progress scales as 0(V^), 
which matches the prior results in the literature. Compared 
with the spatial-reuse multi-hop protocol in 1 1 1 recently proposed 
by Baccelli et ah, in our new protocol the amount of nodes 
involved and the calculation complexity for each relay selection 
are reduced significantly, which is attractive for energy-limited 
wireless ad hoc networks (e.g., wireless sensor networks). 

I. Introduction 

In wireless ad hoc networks, each node may serve as the 
data source, destination, or relay at different time instants, 
which leads to a self-organized network. Such a decentralized 
structure makes the traditional network analysis methodology 
used in centralized wireless networks inadequate. In addition, 
it is hard to define and quantify the capacity of large wireless 
ad hoc networks. In the seminal work |2|, Gupta and Kumar 
proved that the transport capacity for wireless ad hoc networks, 
defined as the bit-meters pumped every second over a unit 
area, scales as &(^/n) in an arbitrary network, where n is 
node density. In |3|, Weber et al derived the upper and 
lower bounds on the transmission capacity of spread- spectrum 
wireless ad hoc networks, where the transmission capacity 
is defined as the product between the maximum density 
of successful transmissions and the corresponding data rate, 
under a constraint on the outage probability. However, the 
above work only considered single-hop transmissions. 

In [4J, with multi-hop transmissions and assuming all the 
transmissions are over the same transmission range, Sousa and 
Silvester derived the optimum transmission range to maximize 
a capacity metric, called the expected forward progress. Zorzi 
and Pupolin extended Sousa and Silvester's work in |5| to 
consider Rayleigh fading and shadowing. Recently, Baccelli 
et al. [l\ proposed a spatial-reuse based multi-hop routing 
protocol. In their protocol, at each hop, the transmitter selects 
the best relay so as to maximize the eff'ective distance towards 



the destination and thus to maximize the spatial density of 
progress. By assuming each transmitter has a sufficient backlog 
of packets, Weber et al. in (E] proposed longest-edge based 
routing where each transmitter selects a relay that makes the 
transmission edge longest. In [7J, Andrews et al. defined the 
random access transport capacity. By assuming that all hops 
bear the same distance with deterministically placed relays, 
they derived the optimum number of hops and an upper bound 
on the random access transport capacity. 

Most of the above works with multi-hop transmissions (e.g.. 
El, O, and ||7l) assume that each hop traverses the same 
distance, which is not practical when nodes are randomly 
distributed. On the other hand, in [IJ and |6| the authors 
proposed routing protocols with randomly distributed relays; 
but they did not address how to optimize the transmission 
distance at each hop. In this paper, by jointly considering 
the randomly distributed relays and the optimization for the 
hop distance, we propose a selection region based multi-hop 
routing protocol, where the selection region is defined by 
two parameters: a selection angle and a reference distance. 
By maximizing the expected density of progress, we derive 
the upper bound on the optimum reference distance and the 
relationship between the optimum reference distance and the 
optimum selection angle. 

The rest of the paper is organized as follows. The system 
model and the routing protocol are described in Section II. 
The selection region optimization is presented in Section III. 
Numerical results and discussions are given in Section IV. The 
computational complexity is analyzed in Section V. Finally, 
Section VI summarizes our conclusions. 

II. System Model and Routing Protocol 

In this section, we first define the network model, then 
present the selection region based routing protocol. 

A. Network Model 

Assume nodes in the network follow a homogenous Poisson 
Point Process (PPP) with density A, with slotted ALOHA being 
deployed as the medium access control (MAC) protocol. We 
also consider the nodes are mobile, to eliminate the spatial 
correlation, which is also discussed in [IJ. During each time 
slot a node chooses to transmit data with probability p, and 



to receive data with probability I - p. Therefore, at a certain 
time instant, transmitters in the network follow a homogeneous 
PPP (Utx) with density pA, while receivers follow another 
homogenous PPP (YIrx) with density (1 - p)A. Considering 
multi-hop transmissions, at each hop a transmitter tries to find 
a receiver in Urx as the relay. We assume that all transmitters 
use the same transmission power p and the wireless channel 
combines the large-scale path-loss and small-scale Rayleigh 
fading. The normalized channel power gain over distance d is 
given by 

(1) 



G(d)= ^, 
d^ 



where y denotes the small-scale fading, drawn from an expo- 
nential distribution of mean - with probability density function 
(PDF) fy(x) = jd exp(-yux), and a > 2is the path-loss exponent. 
For the transmission from transmitter / to receiver j, it 
is successful if the received signal-to-interference-plus-noise 
ratio (SINR) at receiver j is above a threshold ft. Thus the 
successful transmission probability over this hop with distance 
dij is given by 



Ps = ^r 



pyodr; 
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where / g Ilrx, j e IVrx, I^keUTAU} Pyt^kf'' is the sum 
interference from the simultaneous concurrent transmissions, 
dkj is the distance from interferer k to receiver j, and t] is 
the average power of ambient thermal noise. In the sequel 
we approximate y/ = 0, which is reasonable in interference- 
limited ad hoc networks. From [IJ, the successful transmission 
probability from transmitter / to receiver j is derived as 



Ps = exp [-Aptdij^^ , 



where 
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(3) 



(3a) 



sin(27r/Qf) 
B. Selection Region Based Routing 

Considering a typical multi-hop transmission scenario, 
where a data source (S) sends information to its final des- 
tination (D) that is located far away, and it is impossible to 
complete this operation over a single hop. 

Since we assume that nodes are randomly distributed, relays 
may not be located at an optimum transmission distance as 
derived in |7|. To guarantee a relay existing at a proper posi- 
tion, we propose a selection region based multi-hop routing 
protocol. For each transmitter along the route to the final 
destination, we define a selection region by two parameters: a 
selection angle if and a reference distance Vm, as shown by the 
shaded area in Fig. 1, where the selection region is defined 
as the region that is located within angle ip and outside the 
arc AB with \0B\ = r^. Here, the transmitter is placed in the 
circle center O, IBOC = lAOC = ^/2, and OC points to the 
direction of the final destination. At each hop, the relay is 
selected as the nearest receiver node to the transmitter among 
the nodes in the selection region. 

The reason that we limit the selection region within an angle 
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Fig. 1. Selection region 






Fig. 2. Ideal scenario: position-predetermined relays are equidistantly placed 
along the line from the source to the destination 



(f is explained as follows: In multi-hop routing, a transmission 
is inefficient if the projection of transmission distance on 
the directional line from the transmitter towards the final 
destination is negative, or less efficient if the projection is 
positive but very small. Therefore, here we set a limiting 
angle cp with which each packet traverses at each hop within 
[-^/2,^/2]. 

Compared with the model in |7|, where the authors assume 
that the relays are equidistantly placed along a line from the 
source to the destination as shown in Fig. 2, here our model 
assumes that the intermediate relays are randomly distributed 
over the selection region following a homogenous PPP (Hrx) 
with density (1 - p)A, which is more practical. 

III. Selection Region Optimization 

In this section, we optimize the selection region to derive 
the optimum values for the selection angle (f and the reference 
distance r^ by maximizing the expected density of progress. 

As in d, the density of progress is defined as 

D = pA' Ps ' dcoscp, 

where Ps is the successful transmission probability defined 
in (2), dcoscp is the projection of the transmission distance 
d along the directional line OC. Since the receivers follow 
a homogeneous PPP with density A(l - p), the cumulative 
distribution function (CDF) of the transmission distance d is 
given as 



Fr(d < r) = 1 - exp 



-A(l-p)^(r^-ri) 



.rm< r < oo. (4) 



Since is uniformly distributed over [-(f/2,(f/2], which is 
independent of d, the expected density of progress is given by 

-pAtx^ , 



X COS (pfd(x)d(pdx 
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= ^ap{l-p)^l^,kri\k-''^cxpiA{l-p)'^rlym{£j, 



(5) 



where fd(x) is the PDF of d obtained from (4), k = pt -\- (I - 

oo 

/?)|, t is defined in (3a), and r(|,^r^j = J e'H^'^dt is the 

krl 

incomplete Gamma function. 

To optimize the objective function in (5), we first assume 
that if is constant, and try to derive the optimum value of r^. 
For brevity, in the following derivation we write the objective 
function as E. Setting the derivative with respect to r^ as 0, 
after some calculations we have 

dE I ^ 2\ 

— =exp(i(l-;,)-r„) 



upper bound 

numerical result 



3 A drU,krl) 




= 0, (6) 



where r(|,A:r^j is calculated as 

r 1^, /:r2 J = r 1^ j + V^r„ exp {-krl) - ^erf ( V^r„) . (7) 

Thus we have 

dvi^ kr^\ 

-^ = V^exp {-krl) (1 - Ikri) - V^exp {-kri) 

= -2k''V^cxp{-kri). (8) 

Applying (8) to (6), we obtain 

ri^,krl\A(l-p)^r^- Ik'IV^ exp {-krl) = 0. (9) 

Note that the above is only the necessary condition for 
optimality, given the unknown convexity of objective function. 
However, the global optimum must be among all the roots of 
the above equation, which can be found numerically. Since it 
is difficult to analytically derive the exact solution for r^ from 
(9), we turn to get an upper bound of r^. Since 

T{l,kr'^>\[T{lMl)+T{2Ml)] 

= iexp(-/:r^)(2 + ^r2), 

by (9), we have 

^(1 - p)kipri - A^l^r^ + A{\-p)>f> 0. 
Therefore, 

Ik^l'^ - ■^4k^-2k[A(l-p)<fif 
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Fig. 3. Numerical results and the analytical upper bound for the optimum 
reference distance 



For brevity, let us denote Qxp(Akr^) as e and r(|,^r^j as F. 
With partial derivatives, we have 

^=VAp(1-p) sin (I) r^/2 exp {-Aptrl) 



de dY 
-lAptr^Te + T-— + ^-— 



= 0. 



This holds only if 



^de dY ^ 

Y-— +e-— = lAptTmYe. 



(13) 



(14) 



Since k = pt -\-(l - p)j, there is |^ = |(1 - /?). To simplify 
things, we can then calculate the derivative with respect to k 
instead of (f as 

dE 



— = ^/Ap(l-p)cxp[-Aptrl) 



dk 



. cos(^/2) 3 . /^\._K^ ^ . (^\Lde dY\ 



(10) 



= 0. 
(15) 

Since the factor in e and F related to k and Vm is only kr^, 
thus we get ir„^ = k% and \r^§- = k§. 

(11) Therefore, with (14), we have 

^de dY r^l^de dY \ ^ , ^^ ^^^^ 



^m "^ 



kA(l-p)^ 



(12) 



In Fig. 3, we compare the upper bound of the optimum r^ 
with the numerically computed optimal value when (f = ^n. 
We see that when transmission probability p increases, the 
upper bound becomes tighter. 

Now let's maximize the objective function by jointly opti- 
mizing rm and (f. Rewrite (5) as 

E = ^fAp(l-p)cxp{-Aptri)Yi^,krlyxp{Akrl)k-'^^ sm^^ 



Applying (16) to (15), the following holds: 

cos((fi/2) 3 



r? - 2 sin(|)r'Jr..sin(|)r' W„r. = o. (i7) 



After some calculation, (17) is simplified as 

cot(cp/2)k 3 2 n 

-^-—--^Aptr^ = 0. 



(18) 



Since it is hard to derive close-formed solutions for the 
optimal r^ and (f, respectively, we implicitly use (f and p to 
express the optimal r^ as 
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Fig. 4. The expected density of progress vs. the reference distance and the Fig. 5. The expected density of progress vs. the reference distance and the 
selection angle with p=0.01 selection angle with p=0.05 



3/2-cot(^/2)/:/(l-/7) 



Apt 



(19) 



Note that r^ scales as A~^^^, which intuitively makes sense. 
Since as the density increases, the interferers' relative distance 
to the receiver decreases as VI, it requires a shorter trans- 
mission distance by the same amount to keep the required 
SINR. By applying (19) in (5), we observe that (5) becomes 
N ^/A, where N is 3. constant independent of A. This means that 
the maximum expected density of progress scales as 0( VI), 
which conforms to the results in [2J and lj\. 
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IV. Numerical Results 

In this section, we present some numerical results based on 
the analysis in Section III. We choose the path-loss exponent 
Of as 3, the node density /I as 1, and the outage threshold JS 
as 10 dB. In Figs. 4 and 5, we plot the expected density of 
progress vs. the reference distance r^ and the selection angle 
(f, with p = 0.01 and 0.05, respectively. We see that for each 
p there exists an optimum selection angle and an optimum 
reference range when the respective partial derivatives are zero 
as discussed in Section III. 

In Fig. 6, we plot the optimum selection angle obtained 
numerically vs. the transmission probability. As shown in the 
figure, we see that the increment of transmission probability 
leads to the increase of the optimum selection angle. This 
can be explained as follows: The increment of transmission 
probability means the decrement of the number of nodes that 
can be selected as relays; therefore the selection angle should 
be enlarged to extend the selection region. 

In Fig. 7, we compare the optimum reference distance 
obtained numerically with that derived in (19), where (p is 
chosen optimally as that in Fig. 6. We see that the increment 
of transmission probability leads to the decease of reference 
distance. This can be explained as follows: The increment of 
transmission probability means more simultaneous concurrent 
transmissions such that the interference will be increased; 



Fig. 6. The optimum selection angle vs. the transmission probability p 



therefore the reference transmission distance should be de- 
creased to guarantee the quality of the received signal and the 
probability of successful transmission. 

In Fig. m we compare the performance of our routing 
protocol with that in |[T1, and also with the optimized (with 
optimum angle) nearest neighbor routing shown in Fig. 5 of 
[IJ and the non-optimized (with an arbitrary angle, e.g., |) 
nearest neighbor routing in |8|. 

From Fig. 8, we have the following observations and 
interpretations: l)When p increases, the performance of our 
routing protocol becomes close to that of the nearest neighbor 
routing with an optimum angle. This can be explained from 
Fig. 7 as: When p increases, the reference distance r^ tends 
to be a small value close to zero; thus our routing scheme 
degenerates to the nearest neighbor routing. Furthermore, our 
routing protocol shows much better performance than the non- 
optimized (with non-optimized angle) nearest neighbor rout- 
ing, and this advantage is due to adopting both the optimum 
selection angle and the optimum reference distance. In this 
case, the selection region based routing can also be considered 
as the optimized nearest neighbor routing given the selection 
angle and the reference distance; 2) We see that when at the 
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Fig. 7. The optimum reference distance vs. the transmission probabiHty p 
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Fig. 8. The performance of different routing protocols 



optimum selection angle and the optimum reference distance, 
the optimum transmission probability is approximately 0.05. 
Although in this paper we mainly focus on the optimization of 
the selection region, this observation indicates that there also 
exists an optimum transmission probability with our model, 
which has been discussed in some other prior literature, e.g., 
m, 11, 0, and El. 

V. Discussion on Complexity 

As shown in Fig. 8, for multi-hop ad hoc networks, in terms 
of the performance metric of expected density of process, 
Baccelli et al'^ routing strategy is the best, by design at 
each hop the transmitter chooses a relay that provides the 
maximum value of progress towards the destination. However, 
the computational complexity with this protocol might be high, 
since the transmitter at each hop should compute the successful 
transmission probability Ps\x\ together with the projection of 
transmission distance, and accordingly evaluate the value of 
progress towards the destination for each receiver, further 
choose the one with the greatest value of progress as the relay. 

In our protocol as we see from Fig. 8, when p is small, 
its performance is close to that of Baccelli et al, while the 
computational complexity per hop is reduced significantly: 



1) The nodes involved in the relay selection process are 
limited to a small region. As relays are selected from the 
receiver nodes in the selection region, the number of nodes 
participating in the relay selection is reduced with a ratio iplln 
compared with that in [IJ. 

2) Unlike that in [IJ, where the successful transmission 
probability Ps\x\, the projection of transmission distance, and 
further the value of progress towards the destination for each 
potential relay need to be calculated; we only need to calculate 
the distance between the transmitter and the potential relays. 

Also note that our new protocol could be easily imple- 
mented by deploying directional antennas in the transmitter, 
where the spread angle can be set equal to the optimum se- 
lection angle. In this case not only the network computational 
complexity but also the interference will be reduced, which 
will be addressed in our future work. 

VI. Conclusions 

In this paper, we propose a selection region based multi- 
hop routing protocol for random mobile ad hoc networks, 
where the selection region is defined by two parameters, 
a selection angle and a reference distance. By maximizing 
the expected density of progress, we present some analytical 
results on how to refine the selection region. Compared with 
the previous results in JH, O, and Q, we consider the 
transmission direction at each hop towards the final destination 
to guarantee relay eflftciency. Compared with the protocol in 
im, the optimum selection region defined in this paper limits 
the area in which the relay is being selected, and the routing 
computational complexity at each hop is reduced. 
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